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IITRODUCTIOh 

According  to  optical-depth  calculations  [1,2],  106.7  nm  and  104.8  nm 
argon  resonance  radiation  is  extensively  imprisoned  within  the  argon  induc¬ 
tively  coupled  plasma  (ICP).  This  radiation  trapping  lengthens  the  apparent 

radiative  lifetimes  Tapp  and  would  increase  the  populations  of  the  radiative 
3  1 

P}  and  argon  4s  states  if  the  main  relaxation  pathway  were  indeed 
radiative  rather  than  collisional.  Consequently,  radiation  trapping  must  be 
considered  in  developing  both  a  plasma  model  and  an  analyte  excitation 
mechanism  for  the  ICP. 

As  argon  4s  atoms  decay  radiatively,  the  resonance  radiation  is  re¬ 
peatedly  self-absorbed  and  re-emitted  within  the  argon  stream.  Radiative 
loss  of  the  resonance  energy  occurs  only  as  photons  migrate  to  the  stream 
boundaries.  The  degree  of  imprisonment  can  be  described  in  terms  of  the 
probability  that  a  given  resonance  photon  will  escape  from  the  trap.  This 
probability  is  commonly  called  the  escape  factor  and  is  denoted  g.  The 
system  behaves  as  if  the  radiative  resonance  states  possess  lengthened 
lifetimes.  The  resulting  apparent  radiative  lifetime  tapp  is  longer  than 
the  natural  radiative  lifetime  t  according  to 


T  =  T  /g  . 
app  o 


For  the  argon  resonance  lines  in  the  ICP,  g  values  of  10“^  -  10-^  have  been 
calculated  [1,2]  using  Holstein's  theory  for  a  cylindrical  trap  with  Dop¬ 
pler-broadened  line  profiles  [3,4].  This  degree  of  imprisonment  gives 
apparent  lifetimes  of  approximately  1-10  ms  for  the  argon  ^Pj  and  *P^ 
states.  Experimental  values  for  argon  4s  lifetimes  in  the  ICP,  requiring 
spectroscopic  measurements  at  105  nm,  have  not  yet  been  reported. 
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Radiation  trapping  is  important  in  ICP  modeling  in  several  ways.  Im¬ 
prisoned  resonance  radiation  maintains  a  suprathermal  argon  4s  population. 

(A  rapid  and  complete  collisional  mixing  of  the  quartet  of  4s  states  is 
presumed  [1,2,5].)  In  turn,  the  4s  population  is  of  direct  importance  in 
assessing  the  contribution  of  Penning  ionization  to  the  analyte  excitation 
picture.  Furthermore,  for  plasmas  in  which  the  electron  particle  density  ng 
is  near  the  threshold  for  local  thermal  equilibrium  (LTE),  resonance  radia¬ 
tion  trapping  supports  LTE  by  reducing  the  importance  of  radiative  4s  decay 
relative  to  4s  relaxation  by  electron  collisions  [2,6-8].  Lastly,  if  ap¬ 
parent  lifetimes  are  in  the  millisecond  range  [1],  argon  4s  atoms  would 
permeate  the  argon  plasma  body  and  tailflame  (to  the  extent  that  they  are 
not  collisionally  deactivated). 

Previous  trapping  calculations  assumed  the  resonance  absorption  and 
emission  line  profiles  both  to  be  Gaussian  because  of  Doppler  broadening 
[1,2],  in  which  case  the  escape  factor  depends  inversely  upon  the  ground 
state  particle  density  nQ  [3,4].  In  contrast,  for  strongly  self-absorbed 
lines,  escape  of  radiation  from  a  resonance  trap  would  be  caused  largely  by 
losses  in  the  wings  of  the  line  profile  [3,4].  Thus,  trapping  is  particu¬ 
larly  sensitive  to  the  existence  of  pressure  broadening  and  its  Lorentzian- 
dominated  contributions  to  the  line  profile  wings.  At  nQ  values  where 
pressure  broadening  dominates  the  line  wings,  the  escape  factor  g  reaches  a 
limiting  value.  For  laboratory  conditions  this  high-pressure  limit  is  in 
the  range  10“^  <g<  10“\ 

The  argon  resonance  lines  in  the  ICP  possess  a  significant  pressure¬ 
broadening  component,  casting  some  doubt  on  earlier  calculations  [1,2].  The 
full  width  at  half-maximum  intensity  (FWHM)  for  pressure  broadening  of  the 
104.8  nm  line,  for  instance,  is  calculated  for  the  ICP  to  be  -  3  x  10“^  nm 
[9-11]  compared  with  an  estimated  Doppler  FWHM  of  7  x  10”^  nm  for  a  nominal 
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ICP  gas  temperature  of  5000  °K  [10,12].  Although  this  calculated  pressure- 
broadened  component  is  seemingly  less  than  that  of  the  Doppler  effect,  its 
Lorentzian  nature  influences  the  line  wings  dramatically. 

The  Holstein  model  for  a  cylindrical  radiation  trap,  employed  previous¬ 
ly  for  ICP  modeling  [1,2],  presumes  initial  excitation  of  the  resonance 
states  along  the  axis  of  the  cylinder  [3.4].  In  contrast,  the  actual  exci¬ 
tation  geometry  of  the  ICP  is  annular,  with  4s  excitation  occurring  through¬ 
out  the  plasma  region,  presumably  via  the  complex  combination  of  electron 
(and  atomic)  collisions,  ionization,  recombination,  and  radiative  cascading 
expected  in  a  collision-dominated  plasma  [6]. 

These  circumstances,  and  the  lack  of  direct  experimental  lifetime 
measurements,  encourage  further  development  of  resonance  radiation  trapping 
calculations  for  the  argon  ICP.  In  the  present  study,  such  calculations 
are  provided  and  reveal  that  the  radiation  imprisonment  is  governed  by 
pressure  broadening.  The  resulting  apparent  lifetime  Tapp  is  independent  of 
the  argon  ground-state  particle  density  nQ  and  is  considerably  shorter  than 
earlier  calculations  [1,2]  indicated.  The  importance  of  resonance  radiation 


two  treatments  are  in  essential  agreement,  with  steady-state  escape  factors 
being  20  percent  larger  than  the  Holstein  results  for  intermediate  and  large 
optical  depths  (greater  than  -  20).  These  trapping  models  have  been  ex¬ 
tended  by  Walsh  [15,16]  and  Phelps  [17,18].  A  number  of  experiments  [17-25] 
have  shown  the  agreement  between  experiment  and  theory  to  be  generally 
better  than  15  percent.  Here  we  present  a  brief  review  of  Holstein's  theory 
[3, 4]  as  it  applies  to  the  ICP,  combining  Doppler  and  pressure-broadening 
effects  [16]. 

The  transport  of  radiant  energy  through  an  optically  thick  medium  can 
be  described  [3]  by  the  rate  of  change  of  excited-state  particle  density  n 
at  location  r  according  to 


[r]  / 

t  =  J  Y  n 

trap  volume 


(r')  G(r.r')  dr’  -  yn(r) 


where  Y  is  tq“1  and  G(r,r')dr'  is  the  probability  that  radiation  emitted  at 
location  r'  is  transmitted  to  r  and  captured  (absorbed)  there.  The  factor  y 
n(r')  gives  the  rate  of  emission  from  a  volume  element  dr'  at  r'.  The  first 
term  on  the  right-hand  side  of  equation  (2)  thus  describes  the  rate  of 
formation  of  excited  states  at  r  by  absorption  of  radiation  emanating  from 
all  other  locations  in  the  trap.  The  second  term  simply  describes  the  rate 
of  excited-state  loss  at  r  by  emission. 

If  one  defines  a  transmission  coefficient  T(p)  to  describe  the  proba¬ 
bility  that  a  photon  emitted  at  r'  travels  the  distance  p  =  |r-r'| ,  the 
probability  of  photon  capture  over  a  distance  dp  will  be  proportional  to 
-(3T/3p)dP.  For  a  photon  emitted  isotropically  at  r\  the  probability  of 
that  photon  traveling  within  a  solid  angle  doi  is  du/^-rr.  The  probability  of 
photon  capture  G(r.r')  within  a  volume  element  dr  =  p2dpdu>  located  at  r  is 
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then 


G(r,r*) 


(3) 


Substituting  equation  (3)  into  equation  (2)  provides  the  rate  expression  to 
be  solved.  Non-radiative  creation  or  annihilation  of  excited  states  can  be 
included  in  this  treatment  by  adding  the  appropriate  terms  [15,233. 

Solving  equation  (2)  is  complicated  by  the  fact  that  T(p)  depends  upon 
the  spectral  distribution  profile  of  emission  and  absorption.  T(p)  can  be 
written 

T(p)  *  J  P(v)  exp[-k(v)p]  dv  (4) 

where  the  conventional  transmission  probability  exp[-k(v)p3  is  averaged  over 
the  emission  profile  P(\j);  k(v)  is  the  atomic  absorption  coefficient.  For 
Doppler-broadening,  P  (v)  is  proportional  to  k  (v)  [133.  (Also  see  Appendix 
A  of  [43.) 

Figure  1  illustrates  the  nature  of  T(p).  In  Fig.  1  the  absorption 

coefficient  k(v)  has  been  given  the  Gaussian  form  exp[-(v-v  )23  and  the 

o 

optical  depth  kp*  has  been  assigned  the  value  6  at  the  line  center,  to 
describe  a  Doppler-broadened  line  whose  center  is  strongly  self-absorbed 
(Fig.  lc).  For  the  sake  of  illustration,  the  emission  profile  P(v)  has  been 
set  equal  to  k(v);  for  argon  resonance  lines  in  all  regions  of  the  ICP,  P(v) 
is  calculated  [163  to  be  «  k(v)  and  the  magnitude  of  T(p)  is  accordingly 
smaller  than  suggested  in  Fig.  1.  From  Fig.  1,  the  magnitude  of  T(p) 
clearly  derives  from  the  wings  of  the  line  in  a  manner  analogous  to  that 

•The  optical  depth  kp  is  also  known  as  absorbance  or  optical  density. 


found  in  the  integrated  intensity  profile  for  a  self-reversed  atomic  line. 
In  extreme  cases,  radiation  trapping  can  spread  resonance  lines  over  many 
tens  of  nanometers  [26]. 

Equation  (2)  has  for  solutions  a  family  of  eigenfunctions  of  the  form 

n(r,t)  =  n(r)  exp[-gyt]  (5) 


where  g  is  the  radiation  trapping  factor.  Under  pulsed  excitation  condi¬ 
tions,  the  trap  system  would  rapidly  evolve  to  the  n(r,t)  eigenfunction 
having  the  smallest  gy  eigenvalue.  This  smallest  gy  value  defines  the 
escape  factor  g,  the  apparent  lifetime  Tapp  =  (gy)-1,  and  the  corresponding 
"standing  wave"  distribution  n(r).  The  evolution  of  this  longest-lived 
eigenfunction  is  complete  within  a  period  t  «  (gy)"1.  For  steady-state 
excitation,  this  same  longest-lived  eigenfunction  governs  the  predominate 
standing-wave  distribution,  with  contribution  coefficients  for  higher  eigen¬ 
functions  being  inversely  proportional  to  their  corresponding  eigenvalues. 
Holstein  has  calculated  n(r)  and  g  using  the  Ritz  variational  procedure  [3], 
with  g  taking  the  simple  form 

g  1  80  T(P>  (6) 


to  withir.  3  percent.  The  constant  ^  depends  upon  the  choice  of  absorption 
profile  and  trap  geometry.  In  a  cylindrical  radiation  trap,  gQ  is  1.9  for  a 
Doppler-broadened  profile  and  1.3  for  one  dominated  by  pressure  broadening. 
These  values  are  based  on  Bieberman's  steady-state  solution  [14],  and  are  20 
percent  larger  than  those  of  Holstein  [3,4], 
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For  a  cylindrical  radiation  trap  of  radius  R  in  which  initial  excita¬ 
tion  is  along  the  cylinder  axis,  the  radial  distribution  of  excited  states 
n(r)  is  given  by 


n(r)/n{rso)  =  [1  -*r2/R2] 


(7) 


where  <*  is  approximately  8/9  and  3/4  for  the  Doppler  and  pressure-broadened 
cases  respectively.  Figure  2  shows  n(r)  plotted  versus  r  for  these  two 
cases.  When  the  initial  excitation  is  pulsed,  the  distribution  is  estab¬ 
lished  within  a  period  t  <<  (gy)"1,  i*e*.  t  <<  Tapp. 

Walsh  [16]  has  combined  the  effects  of  Doppler  and  pressure  broadening 
under  an  approximate  absorption  coefficient  k(v): 

k(v)  a  kQ[exp(-x2)  ♦  ax2/ir^]  (8) 


where 


X  »  (V  -  Vo)  Xo/v 


o 


a  =  (Y  +  Yc)  Ao/4irvo 
ko  =  Ao  gzYn/Sir*  glvo 

Yc  "  4e2  ™ono/3mc 


1  /2 

vQ  is  the  average  particle  velocity  (2RT/M)  ,  ^  is  the  central  wavelength 

of  the  transition,  f  is  the  oscillator  strength,  and  g2  and  g1  are  the 
degeneracies  of  the  excited  and  ground  states  respectively;  yc  is  the 
average  collision  rate  for  the  excited  state,  kQ  is  the  absorption  coef- 


ficient  at  the  line  center,  nQ  is  the  number  density  of  ground-state  (ab¬ 
sorbing)  species,  m  is  their  mass,  and  e  is  the  unit  charge. 


Using  equations  (4)  and  (8),  T(p)  becomes 
T(P)  »  Td  exp(-,rTCDV4Tc’)  . 


where 

TD  =  [kop(7rlnkoP)l/2jl  ’  Tc  =  [a/TrV^^p] 1/2 


(9) 


ana 

TCD  =  2a/7r(ln  V)V*  • 

TD  and  Tc  are  the  transmission  coefficients  for  pure  Doppler  and  pressure 
broadening,  respectively,  and  TCD  is  the  coefficient  for  pressure-broadened 
emission  and  Doppler-broadened  absorption  profiles.  Eg  is  the  error  inte¬ 
gral  [27]  and  is  1.0  for  the  argument  values  1^2Tqd/2Tc  encountered  in  our 
calculation.  Equation  (9)  reduces  to  Tg  for  large  values  of  the  damping 
coefficient  a_  and  to  Tq  as  a  vanishes,  and  agrees  with  experiment  to  within 
10  percent  [16].  The  optical  depth  kQp  pertains  to  the  line  center.  For  a 
cylindrical  trap  with  axial  excitation,  p  becomes  the  cylinder  radius  R. 

RESULTS 

Cylindrical  trap-axial  excitation.  The  transmission  coefficient  as  a 
function  of  trap  radius,  T(R),  was  calculated  using  equation  (9)  for  a 
cylindrical  trap  of  radius  R  =  0.9  cm  containing  argon  at  1  atmosphere 
pressure  and  over  a  300-10,000  °K  range  of  gas  temperature.  Axial  excita¬ 
tion  was  assumed.  The  temperature  dependence  of  T(R)  derives  simply  from 
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the  variation  in  particle  velocity  (vQ)  and  ground-state  number  density 
(nQ).  Under  these  conditions,  the  optical  depth  kQR  always  exceeds  10-\ 

The  transmission  coefficient  T(R)  was  calculated  to  be  6.3  x  10-1*  for  all 
temperatures  and  for  both  argon  resonance  lines,  despite  their  differing 
oscillator  strengths  (f  is  0.07  and  0.28  for  the  106.7  and  104.8  nm  lines, 
respectively  [9].)  This  calculated  T(R)  value  is  the  high-pressure  limit  Tc 
(cf.  Eq.  9).  For  nominal  ICP  gas  temperatures,  the  first  (Doppler-domi¬ 
nated)  term  of  equation  (9)  contributes  significantly  only  as  the  number 
density  nQ  falls  below  lO1^  cm“3,  which  is  a  full  order  of  magnitude  below 
nQ  in  the  atmospheric  ICP. 

Under  these  conditions  and  using  equation  (6),  the  escape  factor  g  is 

calculated  to  be  8.2  x  10-1*,  giving  Tapp  values  of  10  V13  aid  2.4  us  for  the 
1  1 

P-j  and  P-|  states,  respectively.  These  values  agree  closely  with  direct 
lifetime  measurements  for  pulsed  electron-impact  argon  excitation  at  300  °K 
[28].  The  uncertainty  in  the  trapping  calculation  is  estimated  to  be  15 
percent. 

Cylindrical  trap-ICP  conditions.  Unfortunately,  the  ICP  differs  from 
the  foregoing  Holstein-based  cylindrical  trap  model,  with  its  axial  excita¬ 
tion,  because  of  the  annular  argon  4s  excitation  in  the  ICP  and  the  presence 
of  electron-collision-induced  4s  excitation  and  relaxation. 

Walsh  [15]  addresses  both  of  these  factors  for  the  case  where  electron- 
impact  excitation  permeates  the  trap,  causing  n(r)  to  be  constant  over  the 
trap  radius  rather  than  following  equation  (7)  (cf.  Fig.  2).  By  Walsh's 
calculation,  the  escape  factor  is  increased  by  19?  when  an  atom  experiences 
many  impact  excitations  and  deactivations  during  the  photon  imprisonment 
period  (Tapp).  This  model  rests  upon  near  thermal  equilibrium  between 
radiative  and  impact  processes,  with  a  steady  state  of  excitation-deacti va- 


tion.  This  situation  pertains  more  closely  to  the  ICP  annular  region  than 
to  the  decaying  plasma  of  the  central  channel  and  tail  flame.  In  these 
latter  regions,  xapp  would  be  expected  to  decrease  substantially  as  elec¬ 
tron-impact  deactivation  leads  to  the  loss  of  energy. 

From  Walsh's  calculation,  Tapp  for  the  ICP  annular  region  is  probably 
best  taken  as  8.0  (^P^  and  1.9  ys  (^P.,). 

There  is  ample  evidence  that  the  4s  states  are  equilibrated  within  a 
period  near  10~^s  in  a  collision-dominated  argon  plasma  such  a3  the  ICP 
[2,5].  It  is  then  plausible  to  assume  for  plasma  modeling  that  the  4s 
quartet  shares  a  common  apparent  lifetime  for  the  radiative  loss  of  reso¬ 
nance  energy.  Because  the  3p^  and  ^  states  have  equal  degeneracies,  their 

radiative  decay  rates  are  simply  additive  (t-1  ,,  =  ti-1  +  t?-1)*  Thus 

overall 

the  4s  apparent  radiative  lifetime  effectively  becomes  1.6  ys  for  the  annu¬ 
lar  region,  its  enclosed  central  channel,  and  the  tail  flame,  all  of  which 
share  the  cylindrical  trap  geometry  with  off-axis  4s  excitation. 

Effect  of  argon  eximers  on  radiation  trapping.  The  ICP  is  sufficiently 
dense  in  argon  to  consider  the  formation  of  Ar2*  eximers  from  argon  4s 
atoms,  with  corresponding  eximer  emission  at  126  nm  and  in  the  first  con¬ 
tinuum  to  the  red  of  the  106.7  nm  resonance  line  [28,30].  This  eximer 
emission,  which  cannot  be  absorbed  by  ground-state  argon  atoms,  would  be 
lost  from  the  radiation  trap  and  would  decrease  the  apparent  lifetime  Tapp 
accordingly  [28].  Importantly,  the  ICP  background  spectrum  shows  no  emis¬ 
sion  at  126  nm  [31  ]. 

The  rate  of  Ar2*  formation  in  the  ICP  can  be  estimated,  taking  1.6  x 

op  6—1  i 

10--3*  cm°  s"  as  the  overall  three-body  rate  constant  for  forming  Ar2  from 
Ar  (^P^,  3p.|,  ^Pg)  300  °K  [32].  Adjusting  the  rate  constant  for  a  gas 
temperature  of  5000  °K  [33]  and  converting  to  an  equivalent  pseudo-first- 
order  constant  gives  a  value  of  3  x  101*  s"1  for  an  argon  number  density  nQ  = 
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1.5  x  101*  cm”^.  This  rate  is  an  order  of  magnitude  slower  than  T-i  . 

app 

However,  as  the  gas  temperature  drops  below  about  2000  °K,  the  calculated 
exitner  formation  rate  surpasses  T~pp*  Consequently,  Ar*2  is  most  likely 
present  in  the  ICP  argon  coolant  sheath  and  might  affect  radiation  trapping 
there,  but  probably  does  not  appear  to  a  significant  degree  in  the  central 
plasma  region  or  in  the  sample-containing  area.  Trap  loss  from  eximer 
emission  is  thus  judged  insufficient  to  shorten  appreciably. 

DISCUSSION 

The  degree  to  which  resonance-radiation  imprisonment  maintains  the 
argon  ICP  depends  primarily  upon  the  electron  number  density  ne  and  on  the 
relative  importance  of  electron-collisional  and  radiative  processes.  The 
relative  importance  of  these  factors  can  be  evaluated  with  the  classic 
Bates,  Kingston  and  McWhirter  (BKM)  col lisional-radiative  model  of  a  recom¬ 
bining  plasma  [6].  This  model  is  a  detailed  accounting  of  the  various 
collisional  and  radiative  processes  which  contribute  to  the  recombination  of 
electrons  and  ions  in  a  dense  plasma,  with  subsequent  cascading  of  excited 
states  to  the  atomic  ground  state.  The  overall  net  rate  of  recombination, 
or  rate  of  decay  of  the  plasma  state,  is  expressed  in  terms  of  a  simple  two- 
body  ion-electron  recombination  process,  incorporating  the  consequences  of 
the  complex  set  of  competing  collisional  and  radiative  excitation  and  re¬ 
laxation  events.  The  ICP  would  be  expected  to  exhibit  "recombining"  behav¬ 
ior  [6]  in  the  region  downstream  from  the  load  coil,  but  not  necessarily  in 
the  central  analyte-containing  zone. 

In  the  BKM  model,  net  recombination  is  slowed  by  resonance  radiation 
trapping  to  the  extent  that  the  elevated  resonance  state  population  leads  to 
increased  reionization.  The  effect  of  trapping  is  significant  if  ne  is 


below  a  threshold  value  that  depends  weakly  upon  the  electron  temperature 
Te.  Above  this  threshold,  ne  is  sufficiently  high  that  collisional  pro¬ 
cesses  dominate  excitation  and  relaxation,  and  the  radiative  role  is  negli¬ 
gible.  For  argon  plasmas  at  atmospheric  pressure,  the  nfi  threshold  appears 
to  be  -  2  x  1016  cm"3,  based  upon  the  condition  for  local  thermodynamic 

equilibrium  (LTE)  in  arc  studies  [2,34]  in  which  Tg  was  7000-9000  °K.  Mea¬ 
surements  of  ne  in  the  ICP  generally  give  values  in  the  range  5  x  10^  to 
1  x  10^  [35-37].  Consequently,  one  can  assume  that  recombination  in  the 
ICP  is  retarded  by  resonance  radiation  trapping,  with  the  effect  being 
greatest  at  the  plasma  outer  surface,  in  the  sample  stream,  and  in  the  tail 
flame,  where  ne  is  lowest.  BKM  calculations  for  a  hydrogen  plasma  [6]  show 
that,  as  ne  drops  a  decade  below  the  LTE  threshold  value  at  Te  -  4000  °K, 

resonance  radiation  trapping  slows  recombination  by  only  10  percent.  There¬ 
fore,  the  retardation  in  recombination  should  not  be  dramatic  except  perhaps 
in  the  very  uppermost  region  of  the  plasma  tail  flame. 

The  effect  of  resonance  radiation  imprisonment  upon  the  argon  4s  popu¬ 
lation  is  illustrated  by  the  calculation  of  Giannari3  and  Incropera  for  an 
atmospheric  argon  arc  plasma  [7],  also  performed  using  the  BKM  recombination 
model.  For  10^  £  ne  10  ^  cm"  3  and  7500  _<  Te  <  9500  °K,  complete  reso¬ 
nance  radiation  trapping  causes  the  4s  states  to  be  overpopulated  by  10 
percent  compared  to  the  population  predicted  by  the  Saha  equilibrium  of  Ar+ 
with  the  assumed  number  of  free  electrons.  Without  trapping,  the  4s  states 
are  calculated  to  be  underpopulated  by  50  percent,  compared  to  the  Saha 
value,  as  a  consequence  of  radiational  disequilibrium  induced  by  the  large 
4s  -  3p  transition  probability.  For  ne  the  2  x  lO1^  cm"3  LTE  threshold 
described  above,  the  BKM  calculation  gives  Saha-based  4s  populations, 
whether  or  not  the  resonance  radiation  is  imprisoned. 

In  the  ICP,  the  populations  of  the  4p  (and  higher)  states  should  not  be 
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significantly  affected  by  resonance  radiation  trapping,  except  perhaps  in 
regions  of  lowest  ng.  Even  in  electron-lean  regions  of  the  ICP,  the  4p  and 
higher  states  would  be  expected  to  approach  equilibrium  with  the  4s  states. 
If,  in  these  regions,  the  4s  population  were  elevated  by  radiation  trapping, 
higher-energy  states  would  be  similarly  affected. 

For  axial  4s  excitation  and  in  an  A r  environment  dominated  by  radiation 
trapping,  the  4s  radial  population  distribution  n4a(r)  would  follow  curve  b 
of  Fig.  2.  However,  from  the  BKM  model  [7],  n^Cr)  in  the  ICP  should  be 
determined  primarily  by  electron  collisions  in  regions  of  large  ne.  Ac¬ 
cordingly,  in  the  annular  skin  region  of  the  ICP,  where  ng  is  high,  n^Cr) 
should  resemble  the  radial  ne  distribution  (with  allowance  for  radial  varia¬ 
tions  in  electron  energy  distribution)  and  would  likely  be  within  10  percent 
or  so  [7]  of  the  LTE  4s  population  dictated  by  Te  and  ng.  However,  radia¬ 
tion  imprisonment  in  regions  of  low  nft  (e.g.  in  the  aerosol  channel)  would 
have  a  significant  impact.  In  these  regions  nna(r)  rises  above  that  pre¬ 
dicted  by  LTE. 

In  this  regard,  argon  carrier  gas  in  the  central  aerosol  channel  acts 
as  a  resonance  radiation  sink  as  it  flows  upward  into  the  plasma  core 
region.  In  the  aerosol  channel,  n^  should  reach  the  steady-state  value 
prescribed  by  radiation  imprisonment  within  a  period  t  «  xapp  [33.  For 
Tapp  =  1*6  we  ta^e  this  equilibrating  period  to  be  <  10“^s.  For  a 
stream  speed  of  100  m/s  near  the  sample-jet  orifice  [38],  imprisonment 
equilibrium  is  therefore  established  immediately  as  the  carrier  argon  enters 
the  plasma  region.  The  radiative  transfer  which  establishes  this  equilib¬ 
rium  and  n^g  which  results  from  it  are  independent  of  Tg. 

Radiation  trapping  clearly  predicts  that,  in  the  central  channel,  n 


exceeds  the  value  obtained  from  the  central  channel  n.  and  T_  under  LTE. 


The  absolute  magnitude  of  the  central  channel  4s  population  cannot  be  accu¬ 
rately  determined  from  the  trapping  theory.  However,  It  will  certainly  not 
exceed  the  highest  annular  value;  an  annular  excitation  temperature  of  10000 
-  12000  °K  gives  4  x  10^3  £  4  x  10^  cm“3  when  the  4s  quartet  is 

considered  a  single  level  having  12-fold  degeneracy  and  an  energy  of  11.7 
eV.  In  an  argon  arc  with  an  excitation  temperature  comparable  to  the  ICP, 
Shindo  and  Imazu  1 39 3  found  that,  below  ne  -  2  x  10^  cm“3f  is  indepen¬ 
dent  of  ne,  largely  because  of  radiation  trapping.  Their  measurements 
indicate  that  the  total  4s  population  density  is  approximately  1  x  1 0 1  ^  cm” 3 
in  a  non-LTE  argon  plasma  with  trapped  resonance  radiation.  Uchida  and 
coworkers  [37]  have  recently  reported  relative  radial  ^P2  distributions  in 
the  ICP  which  indicate  that  the  central  channel  n4s  is  substantially  smaller 
than  in  the  annular  region,  with  the  radial  gradient  diminishing  with  eleva¬ 
tion  and  disappearing  above  an  elevation  of  approximately  15  ran  from  the 
load  coil. 


CONCLUSIONS 

Imprisonment  of  argon  resonance  radiation  in  the  ICP  is  limited  by 
losses  caused  by  pressure  broadening  of  the  absorption-emission  line  pro¬ 
files.  The  theoretical  radiation  escape  factor  g  calculated  for  a  cylindri¬ 
cal  ICP  radiation  trap  with  annular  argon  4s  excitation  (caused  predominate¬ 
ly  by  electron  collisions),  has  a  lower  limit  of  7.9  x  10"1*.  Apparent 
radiative  lifetimes  Tapp  are  found  to  be  8  us  and  1.9  us  for  the  3p1  and 
1P1  argon  resonanoe  states,  respectively.  The  quartet  of  4s  states,  which 
are  assumed  to  be  completely  interconverted  by  electron  collisions  under  ICP 
conditions,  share  an  overall  apparent  radiative  lifetime  Tapp  ^  us. 


This  apparent  lifetime  is  independent  of  the  argon  gas  temperature  and 
relatively  insensitive  to  whether  the  plasma  state  is  dominated  by  radiative 
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or  collision-induced  processes. 

Radiation  trapping  increases  the  argon  4s  atom  density  by  an  esti¬ 
mated  10  percent  over  that  predicted  by  LTE  in  the  collision-dominated 
regions  of  the  plasma  wherever  ne  >  l  x  1015  cm-3  (e.g.  in  the  plasma 
"fireball").  Substantially  larger  deviations  from  LTE  are  predicted  in 
regions  (e.g.  the  central  channel)  where  nfi  falls  well  below  this  value. 
Radiation  trapping  is  predicted  to  retard  electron-ion  recombination  by  10 
percent,  which  contributes  to  the  ambient  ng  and  argon— ion  concentrations. 
The  trapping  supports  a  suprathermal  njjg  in  the  central  aerosol  channel  at 
all  elevations  in  the  plasma. 

From  these  considerations,  radiation  trapping  would  seem  not  to  be 
responsible  for  sustaining  the  plasma  well  beyond  the  load-coil  region,  as 
earlier  postulated  by  Blades  and  Hieftje  [1].  The  trapping  lifetime  of  1.6 
us  reported  here  would  support  the  discharge  no  more  than  0.05  mm  beyond  the 
energy-addition  zone,  assuming  a  gas  velocity  of  30  m/s.  Clearly,  another 
mechanism  for  ICP  maintenance  must  be  sought;  plasma  modulation  experiments, 
now  underway  in  our  laboratory,  are  directed  at  identifying  the  responsible 
process . 
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FIGURE  CAPTIONS 

Figure  1.  Illustration  of  the  nature  of  the  transmission  coefficient  T(p) 
for  a  strongly  self-absorbed  Gaussian  line:  a.  Gaussian  profile 
for  absorption  coefficient  k(v)  and  emission  profile  P(v).  b. 
Transmittance  exp[-k(v)p].  c.  P(v )  exp[-k(v)p].  T(p)  is  the 
area  under  curve  c  according  to  equation  (4).  For  b,  k  is  taken 
as  6  at  the  line  center  to  give  strong  self-absorption.  Note 
that  for  argon  resonance  radiation  in  the  ICP,  P(v)  <<  k(v), 
and  T(p)  is  appropriately  smaller  than  portrayed  here. 

Figure  2.  Calculated  steady-state  radial  distribution  of  excited  states 
n(r)  in  a  cylindrical  radiation  trap  of  radius  R  with  initial 
excitation  along  the  cylinder  axis  a.  Doppler  line  profile,  b. 
Pressure-broadened  profile.  Both  distribution  curves  have  been 
normalized  to  the  same  peak  value. 
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